Cold and Warm Denaturation of Hydrophobic Polymers. 
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We introduce a polymer model where the transition from swollen to compact configurations is due 
to interactions between the monomers and the solvent. These interactions are the origin of the 
effective attractive interactions between hydrophobic amminoacids in proteins. We find that in the 
low and high temperature phases polymers are swollen, and there is an intermediate phase where 
the most favorable configurations are compact. We argue that such a model captures in a single 
framework both the cold and the warm denaturation experimentally detected for proteins. Some 
consequences for protein folding are discussed. 
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Modeling polymers and polymer collapse, beyond be- 
ing a challenge of great theoretical interest, is extremely 
important for many different applications jj| . In particu- 
lar, in connection to the protein folding problem, the col- 
lapse of polymers has gained a special status in statistical 
physics. The simplest models of proteins describe them 
as polymers whose building monomers (amminoacids) 
can be of two different types, either hydrophobic (H) or 
polar (P) (such a simplification was introduced in the so- 
called HP model Q). The protein folds through attrac- 
tive H-H interactions. (Homo)Polymers with attractive 
monomer-monomer interactions are known to undergo a 
phase transition from a high-temperature swollen phase 
to a low-temperature compact one. The presence of at- 
tractive interactions between hydrophobic amminoacids 
justifies then the collapse of the protein below the folding 
temperature. In particular, trying to maximize the num- 
ber of H-H contacts, most hydrophobic monomers will 
be found in the core of the protein, in agreement with 
experimental observations |3| . 

It is important to understand the origin of the attrac- 
tive H-H interaction: as their name suggests, hydropho- 
bic amminoacids do not like being in contact with water. 
As a consequence, hydrophobic amminoacids in solution 
tend to aggregate in order to minimize the area exposed 
to water. This tendency to aggregate can be modeled as 
an attractive interaction. In proteins, hydrophobic am- 
minoacids try to hide away from water, burying them- 
selves in the core of the protein (whose surface is mainly 
made of polar, hydrophilic amminoacids). 

At temperatures much lower than the folding transi- 
tion, the HP model (and models derived from it) predicts 
that proteins stay in a compact state, with a greater and 
greater probability to be in their ground state. This prin- 
ciple has also been used to classify different amminoacid 
sequences as good or bad folders (from a thermodynam- 
ical point of view) according to e.g. the gap between the 
ground state and the first excited states, or the unique- 
ness of the ground state on a compact configuration Q . 

Recently, nonetheless, there has been a growing evi- 



dence for the so-called cold denaturation of proteins: at 
low temperatures proteins such as /3-Lactoglobulin A, 
myoglobulin, apomyoglobulin, Ribonuclease A, and Es- 
cherichia Coli's HPr Q unfold to a swollen configura- 
tion fi. Cold denaturation is assumed to be a general 
property of all globular proteins . These experimental 
findings are clearly incompatible with the predictions of 
the HP model. 

In this Letter we propose a mechanism driving the 
collapse of a polymer, mimicking the nature of the hy- 
drophobic effect. Furthermore, we find within a single 
framework both a cold and a warm collapse transitions, in 
agreement with experimental observations for hydropho- 
bic amminoacids. The model is built on the present un- 
derstanding of the microscopic organization of the wa- 
ter/hydrophobic molecules system. 

Hydrophobic molecules are essentially non-polar enti- 
ties, weakly, if at all, interacting with water molecules. 
Therefore, there is no true repulsion between water and 
hydrophobic molecules. We can get an insight in the 
energetics of water/hydrophobic molecules systems from 
experiments. In the case of pentane M (hydrocarbons are 
in general hydrophobic), one finds an enthalpy of trans- 
fer AH a ^ s from the aggregate phase to the aqueous so- 
lution phase that is strongly temperature dependent: at 
low temperatures AH a ^ s < 0, and at high temperatures 
AH a ^ s > 0, showing that at low temperatures the solute 
phase is energetically more favourable than the aggregate 
one; at high temperatures the situation reverses. 

Microscopically, a hydrophobic molecule in solution 
displaces water molecules that, at a low enough tem- 
perature, build an ice-like cage around the hydrophobic 
molecule, giving origin to a structure that is energetically 
more favorable than bulk, liquid water || . This energetic 
gain gives the major contribution to the minimization of 
the free energy, and as a consequence the system tries 
to maximize the number of cages by exposing as many 
hydrophobic molecules as possible, hence the cold denat- 
uration of proteins and AH a ^ s < 0. At higher temper- 
atures water around hydrophobic molecules can no more 
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form the cages, but due to the steric constraints imposed 
by the molecule, they cannot even fully exploit the ener- 
getically favorable hydrogen bonding of bulk, liquid wa- 
ter. As a result, disordered water molecules around hy- 
drophobic amminoacids are energetically less favorable 
than bulk water (AH a ^ s > 0): proteins try therefore to 
hide their hydrophobic parts in their core. This is what 
is usually referred to as the hydrophobic effect, and it is 
the consequence of a complex interplay between energy 
minimization and entropy maximization. Unfortunately, 
the detailed behavior of water at low temperatures is 
still far from being well understood from a microscopic 
point of view; the nature of the hydrogen bond itself is 
still a matter of controversy and of deep investigations 
|1. Therefore it is still not possible to give a detailed 
description of the hydrophobic amminoacid-water orga- 
nization and of the energies associated to the cages and 
to the disordered configurations. 

We model the polymer as a self-avoiding walk (SAW) 
on a lattice. In every lattice site (except those occupied 
by the polymer) there is a Potts variable with q states 
(labelled for convenience from to q — 1), representing a 
group of water molecules in q different collective states. 
We associate the state q = to the cage configuration, 
energetically favorable when water is in contact with the 
polymer, and the remaining q — 1 states to disordered, 
unfavorable, configurations. The Hamiltonian of the sys- 
tem is 

N 

H = E E '(- J ^.o + *(i - *»„o)) • (i) 

z=l jn.n.i 

The first sum runs over the N monomers of the poly- 
mer. The second primed sum runs over the water occu- 
pied nearest neighbors of each monomer. The interaction 
constants J and K (both positive) represent respectively 
the energies of the cage configurations and of the dis- 
ordered ones with respect to bulk water. There are no 
monomer-monomer interactions. 

Some kind of water-protein interaction with a similar 
description of the water degrees of freedom was already 
proposed in |L0| with the same motivations. Yet, in that 
model, the water-protein interaction was introduced in a 
"mean field" fashion, and moreover protein folding was 
described by collective hierarchical variables hiding the 
microscopic description, that we believe fundamental to 
give qualitative, but also quantitative, predictions. 

Starting from (|l|) we can write the partition function 
of the system as 

Z n =J2 z n(C) (2) 

c 

where Zn(C) is the partition function associated to a 
single configuration C. It is important to observe that 
the maximum number of water sites in contact with the 



polymer is M — 2(d — 1)N + 2 for a hypercubic lattice 
in d dimensions. For general polymer configurations the 
number of contacts is smaller than M. Nonetheless, all of 
these M water sites must be taken into account in order 
to give the correct weight to all the Zn(C)'s. Moreover, 
the way Hamiltonian (|IJ) has been written implies that 
a single water site can be counted more than once. It is 
then possible to write a fairly simple expression for the 
configuration partition function: 

z 

Z N (C) = q n °W[[Y l (C) (3) 

i=i 

with 

Y l (C) = ((q-l)e-f ,lK + e^ J ) n ' (C) (4) 

where n; (C) is the number of water sites with I polymer 
contacts and z is the coordination number of the lattice. 
As usual, P = 1/ksT and we take Icb — 1. Analogously, 
it is possible to write similar expressions (even though a 
little more complicated) for the internal energy U (C) of 
a configuration (and then calculate the internal energy 

We analyze the thermodynamic behavior of the system 
in d = 2 by means of exact enumeration techniques on the 
Manhattan lattice for polymers of length up to N = 25 
monomers. The Manhattan lattice is a two-dimensional 
lattice on which rows (columns) are alternately left /right 
(up/down) oriented. The physics of polymer collapse in 
the well-known case of attractive monomer-monomer in- 
teractions is the same on the Manhattan lattice as on 
regular two-dimensional lattices (there is still some de- 
bate on non-universalities in the values of some exponents 
UPl)- We are therefore confident that our exact enu- 
meration on the Manhattan lattice describes the correct 
behavior of the polymer on regular Euclidean lattices. 

We calculated the specific heat per monomer of the sys- 
tem as the derivative of U with respect to the tempera- 
ture T, see Fig|| We used K/J — 2 and q — 16000 (both 
K and the temperature can be normalized with respect to 
J). In Fig. [j] two peaks of the specific heat appear, corre- 
sponding to temperatures Tq and 7\y. By direct inspec- 
tion of the values of the configuration partition function 
in eq.(^), we find that below Tc the most probable con- 
figurations are swollen, maximizing the number of water- 
polymer contacts. Between Tc and Tw the polymer 
collapses, and the most probable configurations are the 
compact ones (in the particular cases of N — 16, 25 we 
verified that they correspond to square configurations). 
Finally, for T > T\y the weight of compact and swollen 
configurations becomes comparable and the warm denat- 
uration takes place. 

A number of questions arise for this problem: as it can 
be seen from Fig.|], the low-temperature peak in the spe- 
cific heat is stable against the length of the polymer, a 
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clear indication of a true phase transition. On the con- 
trary, the high-temperature specific heat peak seems to 
flatten as the length of the polymer is increased, indicat- 
ing that the high temperature denaturation transition 
could disappear as the length of the polymer increases. 
This scenario needs further investigations. Moreover, the 
behavior of the two peaks as a function of K/J is such 
that they coalesce when K/J = 1. The height of the 
warm denaturation peak depends also on q, increasing as 
q increases, as shown in Fig.^ (as a consequence, the flat- 
tening of the warm denaturation peak depends on the 
interplay between q and N). Also the cold and warm 
denaturation temperatures depend on q. In the inset 
in Fig.|| we temptatively show that the dependance of 
TV on q follows a power-law: the best fit, quite stable 
adding points on both extrema (q small and large) is 
T w ~ 0.22 + lA7q-° 17 . 

Hydrophobic amminoacids are modeled through an ef- 
fective mutual attractive interaction. In this way the 
solvent degrees of freedom can be neglected. Here we 
show that the Hamiltonian in eq.(0) and the correspond- 
ing configuration partition function in eq.(^) give rise 
to an effective monomer-monomer attractive interaction. 
Let us consider the two polymer configurations C and 
C represented in Fig.0. Configuration C is character- 
ized by having a monomer-monomer contact, as opposed 
to configuration C. If we forget about the solvent, and 
we attribute the difference of the two partition functions 
Z(C) and Z(C) to a monomer-monomer interaction en- 
ergy e, we can write 



1 Z(C>) 



(5) 



The energy e is represented in Fig.|] as a function of the 
temperature. The low and high temperature limits can 
be easily computed as 



2J 
—2\K — 



J + K 

q 



t -> o 

T — ► oo 



(6) 



The effective monomer-monomer interaction energy is 
therefore repulsive for small temperatures, and attrac- 
tive for high temperatures, mimicking an effective attrac- 
tive H-H interaction. Yet, through this model, it is ex- 
tremely clear that the hydrophobic attractive interaction 
is present only when hydrophobic monomers are in con- 
tact with water. We can picture a situation where only a 
few monomers are hydrophobic. In this case they try to 
hide in the core of some globule. Within the core there 
are no more interactions between them, if not residual 
ones such as Van der Waals forces. As a consequence, 
in the core of a protein, the hydrophobic interaction can 
not justify the maximization of the number of contacts, 
neither it can discriminate between the native state of 



the protein and the molten globule state, where the hy- 
drophobic protein core is believed to be in an amorphous 
state. 

The absence of interactions in the core of the collapsed 
polymer is the cause of the weakness of the warm de- 
naturation transition: for a polymer of N monomers, 
the energy of compact configurations grows like iV 1 / 2 , 
since it is just a perimeter effect. On the contrary, we 
can expect an exponential number of ways to open up 
the compact polymer, and therefore the entropy grows 
like N. On larger and larger length scales the entropic 
term dominates and compact configurations are not sta- 
ble anymore. On the contrary, cold denaturation is a true 
transition: again the entropy grows like N, but in this 
case the energy is extensive too. 

In conclusion, we have introduced a model Hamilto- 
nian for the collapse of a polymer interacting with the sol- 
vent, but without monomer-monomer interactions. Such 
a model mimics (in an extremely simplified way) the be- 
havior of hydrophobic polymers in water: as a result we 
find the presence of both a cold and a warm denatura- 
tion temperatures, as experimentally found for proteins 
(whose folding is believed to be induced by hydropho- 
bic interactions). Between Tyy and Tc the polymer col- 
lapses to compact configurations. Indeed the model tries 
to capture some of the essential physics of the water- 
hydrophobic amminoacids interaction at a microscopic 
level, and as a byproduct (necessary for consistency), it 
justifies the presence of an effective attractive H-H inter- 
action at high temperatures. Yet, this model rises some 
questions about the reliability of using an attractive H- 
H interaction also in the core of proteins (where water 
is absent) at odds with models presently in use. This 
model further shows that the solvent degrees of freedom 
can, and should, be taken into account in the formu- 
lation of protein models, at least in a simplified way. 
In particular it calls for a much better understanding of 
the physics of water around hydrophobic amminoacids. 
We are presently working on larger computer simulations 
of this model for polymers embedded in lattice different 
from the Manhattan one. We are also considering slight 
modifications of this prototype model to take into ac- 
count a a better description of the water degrees of free- 
dom and of their energetics. Nevertheless, we believe to 
have pointed out the essential ingredients to describe in 
a single framework the processes of cold and warm de- 
naturation. 

We thank P. Bruscolini, L. Casetti, A. Hansen, G. 
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FIG. 1. Specific heat of the system for different polymer 
lengths N = 16, 20, 25. Here K/J = 2 and q = 16000. The 
cold and warm denaturation temperatures, Tc and Tw, are 
marked with arrows. 




FIG. 2. Specific heat vs. temperature for K/J = 4 and dif- 
ferent values of q, namely q = 2™ • 1000 with n — 0, 8, and 
N = 16. In the inset the cold denaturation temperatures (+) 
are shown vs. q (they have already been shifted by the con- 
stant term 0.22); the solid line is the fit JV-0.22 ~ 1.47g~ 017 
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FIG. 3. Two different configurations differing only for a 
monomer- monomer contact. These configurations are used to 
calculate e as in Fig. 3. The water sites that changed their 
number of contacts from C to C' are indicated with their 
number of contacts. 
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FIG. 4. Values of the effective monomer-monomer in- 
teraction energy as computed from (5) and the configura- 
tions in Fig. 2. The ratio K/J — 2 and the values of 
q = 16000, 2000, 1000, 100, 50, 10 have been used. The transi- 
tion between repulsive and attractive energies becomes more 
and more abrupt as q increases. 
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